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Abstract. We prove the existence of solutions for the singularly perturbed 
Schrodinger— Newton system 

with an electric potential V that decays polynomially fast at infinity. The 
solution i/i concentrates, as ?i. — > 0, around (structurally stable) critical points 
of the electric potential. As a particular case, isolated strict extrema of V are 
allowed. 

1. Introduction and statement of the result 
Systems of Schrodinger-Newton type like 

^ ' \ AU + 47r7|?Ap = 

were introduced by R. Penrose in [19 to describe a system in which a mass point 
(of mass m) is placed at the origin, under the effect of the gravitational field, and 
interacts with a matter density given by the square of the wave function -0, which is 
the solution of the Schrodinger equation. In addition, an electric field is generated 
by a potential V. In U is the gravitational potential, 7 — Gm? , G is Newton's 
constant, and h is Planck's constant. 

A few papers have recently appeared in the literature about this (or some related) 
problem, with particular attention to the so-called semiclassical regime h 0. 
In [52], the existence of multibump solutions concentrating (as fi, ^ 0) at local 
minimum points of V is proved, under the assumption that 

inf V(x) > 0. 

This is a broadly-used assumption even in the study of a singular Schrodinger 
equation of the form 

(2) -h'^Au + V{x)u = f{u). 

Many efforts have been made to prove existence results (both for ?i = 1 and ft ^ 0) 
for ([2|) when the potential V is either null in some region (|8j), or tends to zero at 
infinity ([HHH]), or even both things ([S])- In this note, we are interested in the 
second case, under some additional assumptions on the rate of decay at infinity. It 
will be convenient to scale variables in ([1]): 

h ^J^TT^m ' 2m ' 2m 

Then ([1]) becomes 

" ft2AV> - V{x)i: + JJV- = „3 



(3) 



ft2At/ + 47r7|V'p = 
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By classical potential theory, we can solve the second equation of ([3]) with respect to 
U , and plug the result into the first equation. We need to solve the single non-local 
equation 

We set for typographical convenience e = h, u = and W{x) — l/(47r|a;|). We 
will study the singularly perturbed equation 

(4) ~-e^Au + V{x)u = -^{W*\u\'^)u in R^. 

Remark 1 . Formally, the local nonlinearity corresponds to the convolution kernel 
W — S, the Dirac delta distribution. 

We will retain the following assumptions on V: 

(VI) there exist constants Aq, > and there exists a £ [0, 2], such that 



< V{x) < Ai for aU x e 



1 + \x\°' 

(V2) there exists a constant Vi > such that 

sup \yv{x)\ < Vi. 

The second assumption is rather technical, while the first one allows the potential 
V to decay to zero at infinity, although not faster than polynomially. 
A simple change of variables yields 

(5) - Au + V{ex)u= {W 'i'\u\^)u in 

We remark that the homogeneity of W produces cancellation of e in the right-had 
side of ©. 

Remark 2. A similar model has been studied in |16j . with a similar technique. 
However, the equation of [16] does not scale coherently in the non-local term. In 
some sense, the Authors study a system like (O, in which the second equation is 
not singularly perturbed. On the other hand, they allow very general convolution 
kernels W, even without homogeneity, but need a periodic potential V. However, 
we do not know any (non-trivial) explicit convolution kernel that matches all the 
requirements of [16j . in particular the non- degeneracy of the limiting equation (see 
Theorem [S] for a definition of non-degeneracy in our context). 

Remark 3. In the recent paper |12j . the Authors have studied a non-local equation 
like ([4]), under the action of an external magnetic field. The technique is based on 
a compactness analysis of the set of ground-states for the limiting problem. We do 
not believe that this method can be directly applied to the case under consideration. 

We can state our main existence result. 

Theorem 4. Under the assumption (VI) and (V2), suppose that V is smooth 
and that xq € M.^ is an isolated strict local minimum (or maximum) point of V . 
Then for every e small enough, there exists (at least) a solution Ug of ^ such that 
/r3 |Vwep -|- y(ea;)|Mep < oo. This solution is positive and concentrates at xq as 

Remark 5. Our solution is not, in general, a ground-state solution. Even for a single 
Schrodinger equation with decaying potential, ground-state solutions need not exist 
at all (see [2]). However, it is a bound-state solution, since /jja |VMep-|-y(ej:)|Mep < 
oo. 
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Remark 6. It is easy to show that if a bound-state solution concentrates at some 
point xo G M^, then Vy(xo) — 0. This can be done as in [T] or pp] . 

Remark 7. More general existence results can be proved. In particular, we can 
allow V to possess a whole (compact) manifold of non-degenerate critical points. 
We will state this result in the last section of the paper. 

We spend a few words about the technique we employ to solve ([5]). This equation 
has a variational structure, so that we may use tools from Critical Point Theory. 
The elegant approach of [2] does not seem to apply, since it is based on fine in- 
equalities for the local nonlinearity of the equation in order to recover compactness. 
Similar non-local inequalities are not available, as far as we know. The penalization 
method introduced in [71 HH] (a refinement of the celebrated penalization technique 
introduced in jl3p consists of changing the nonlinearity outside some bounded re- 
gion, and again this seems quite hard to adapt, since in our equation the right-hand 
side depends on the values of the unknown itself over the whole space. 

The most flexible approach seems to be based on a perturbation technique in 
Critical Point Theory, introduced and widely applied by Ambrosetti et al. some 
years ago. We refer to the book [3] for the details. In particular, we will modify 
the scheme of [S] and [3] to adapt it to our equation. Loosely speaking, we will 
(1) construct a set of "quasi-solutions" , and (2) construct around this set a finite- 
dimensional natural constraint: our problem will be therefore reduced to the search 
for critical points of a function of three variables. Some difficulties will arise, due 
to the fact that the potential function V can approach zero at infinity. 

Indeed, if V is bounded away from zero, one can work in the usual Sobolev space 
H^{B?). This is what is done in [22] ■ Since V is allowed to decay to zero, the 
integral /j^a V\u\^ need not be finite for every u e H^{M?). It would be natural to 
work in some weighted space, but then the non-local term may be undefined. This 
kind of trouble already arises in single Schrodinger equations, see [2]- We will adapt 
an idea introduced in 0] (see also [5]): a suitable truncation inside the convolution 
will be performed, in such a way that a new functional will be defined smoothly on 
the natural weighted Sobolev space. We stress that this truncation is not used in 
order to confine solutions in any region of the physical space. The concentration 
behavior of our solutions as e — > is not a consequence of this penalization (as it 
is in the scheme of 7 , for example) , but follows from the construction itself. 

Notation. 

(1) D^^^{^^) is the space of all measurable functions w: — > R such that 
/ju3 iVwp dx < Qo. This space is continuously imbedded into L^(R'^), since 
6 is the Sobolev critical exponent in dimension three. 

(2) To save space, we will often write J instead of /jg3 (and also instead of 

(3) The letter C will denote (positive) constants whose precise value is not 
made explicit, and that may vary from line to line. Inside a proof, we may 
also write Ci, C2, C3, . . . , to denote possibly different constants. 

(4) The symbol V and will be freely used to denote first and second deriva- 
tives of functions and functionals, especially when derivatives are thought 
of as vectors (through the Riesz isomorphism). 

(5) Sometimes we will use Landau's symbols: Oe(l) will stand, for example, for 
a quantity that tends to zero as e — > 0. 
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2. The variational setting and background results 
As announced in the Introduction, we will work in the weighted space 

E ^ E, = l^ue D^^'^iR^) \ J V{ex)\u{x)\^ dx < +00^ , 

which is a Hilbert space with respect to the inner product 

{u,v) — {u,v)e ~ / VM-Vu(ia;+ / V{ex)u{x)v{x) dx 

and corresponding norm 

= IIuIIp = / \Vu\'^dx+ / V{ex)\u{x)\'^ dx. 



In the sequel, we will make an extensive use of the following Hardy-Littlewood- 
Sobolev inequality. We refer to |15l I21j for a proof, in a more general setting. 

Proposition 8. Let f e LP{m.^), g G L\M.^), with i + i + ^ = 2. Then 



f{x)9{y) 
\x - y\ 



dx dy 



<c\\f\\M\t 



for some universal constant C > independent of f and g. Moreover, if r, s € 
(1, +00) and < A < 3 satisfy ^ + ^ = 1 + , then 

giy) 



(6) 



-y\ 



■dy 



<C\\g\\ 



For uG E, the double integral /jja^jja W{x — y)u{x)'^u{y)^ dx dy may be infinite. 
To exploit the variational structure of (O, we introduce a smooth function /i € 
C°°(R) such that < ^ < 1, ^(s) = 1 if s < 1, and n{s) = if s > 2. For any 
u G E, we set 

T{x,u) = ^l{c-'\u\■{l+e\x\f), 

where the constants c and 9 will be chosen appropriately later. Then, we define the 
truncation 

Tf;u{x) — T(x, u)u + (1 — T(a;, u)) 



In particular. 



TeU{x) 



u{x), 



{l + e\x\y 



if u{x) < c (1+ e|a;|)~ 



^c(l + £1x1)"" if > 2c(l + e|a;|)"'' 
and we can always choose > so that the functional : i? ^ R defined by 

W{x-y)\TMx)\^\TMy)\^dxdy 



Ie{u)^-\\u\\l 



is of class C^{E). As such, every critical point of such that < 
c(l + e|a;|)~ is a solution of ([S]). 

To proceed further, we need to recall some results about ([5]) with e fixed. 

Theorem 9. There exists a unique radial solution of the problem 

'AU + [/ = (VK * |C/|2) U in M3 



u >0 

U{0) = max2.gR,3 U{x). 



m . 
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The solution U is a strictly decreasing function that decays exponentially fast at 
infinity together with its first derivatives. Moreover, it is non-degenerate in the 
following sense: if <p (z H'^{M.^) solves the linearized equation 

-A(j) + = (M^ * |C/p) (j) + 2{W^ {Ucf))) U, 

then (j) is a linear combination of dU/dxj, j — 1,2,3. 

Proof. The existence of U is proved in [T7]. The second part of the Theorem is 
proved in [H]. □ 

Remark 10. Since the unique solution U of Theorem IH] has a mountain-pass char- 
acterization, it is not hard to refine the statement as foUows: the Euler-Lagrange 
functional associated to the equation solved by U has a second derivative that is 
negative-definite along the direction U, and positive-definite on the orthogonal of 
the subspace defined by the (three) partial derivatives of U itself. 

Remark 11. The exponential decay of the solution U will be used to get high-order 
summability when needed. Many integral estimates would be more involved (and 
even wrong), if U had just a power-like decay at infinity. We refer to the comments 
at the end of the last section. 

For every a > 0, the function Ua{x) — aU{\/ax) is the unique solution of 

^AUa+aUa^ {W *\Ua\')Ua iuM^. 

Given ^ G M^, we set 

Zedx)^Uvieo{^-O^V{ei)u(yvW){x-0) ■ 
Hence z^^^ solves the limiting equation 

(7) -Az,,^ + ^(£0^e,« = (M^*Ne,d')^e,e inM^ 

Although ([5]) cannot be set in H^iJS?) because V can be arbitrarily small, not all 
is lost: Ze^^ is almost a critical point of I^. This is rather standard if V is far from 
zero (see [31 [S]), but it is still true under our assumptions. See Lemma [T^ below. 
Therefore, we will look for a solution of ([5]) near the manifold 

z = z, = {z,,e UeM^ |£?| <i}. 

To this aim, we will perform a finite-dimensional reduction. If P = P^^^ stands 
for the orthogonal projection of E onto {T^^ ^Z)-^, we will first solve the auxiliary 
equation 

(8) PVI,iz,^^+w) =0 

with respect to u; = w{£,^) in some suitable subspace orthogonal to Z. Since we 
can prove that Z = -|- w} is a natural constraint for I^, we will show that 
an isolated local minimum or maximum point xq of V generates a critical point of 
the constrained functional on Z. Thus our main existence result will be proved. 
This perturbation method is a modification of that introduced in [5], and follows 
closely [3]. Since the only difference with [1] lies in the convolution term, we will 
often be sketchy on those arguments that do not involve it, and refer to [3] for the 
details. 
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3. Solving the auxiliary equation 

We recall an analytic estimate on V. We omit the proof, since it appears in [1]. 

Lemma 12. Let a > 0, suppose that V{x) > a|a;|^" for \x\ > 1, and let m > be 
given. Then there exist Eq > and R > such that 

Viex + y)> T-r- 

whenever \x\ > R, e < Eq, y G M'' and \y\ < 1. 

Next, we prove that Z is a manifold of quasi- solutions of We begin with 
a useful fact that will be systematically exploited hereafter. Since z^.j decays ex- 
ponentially fast at infinity, we can (and do) choose c and 9 in such a way that 
TgZe^S, = '^E-i- Therefore, we can simply write z^.j when we work with on Z . 

Lemma 13. There is a constant C > such that ||V/e(zej)|| < Ce whenever 

H\ < 1- 

Proof. Fix any v Cz E. Then 

S/I,{z,^^)v = J {\/z,^^-\/v + V{ex)z,^^v) 

-J W{x - y)zs^^{x)'^Zs^^{y)v{y)dx dy 

and from 
we get 



|VJe(z,,^)i;| < 



[V{ex) - V{eO] z^,^{xfv{x) dx 



We have cleared the convolution term, and the proof reduces now to the arguments 
and the estimates of [U Lemma 3]: we apply the Intermediate Value Theorem to 
the difference V{ex) — V{e^), and we exploit the exponential decay of z^^j to get 
the necessary integrability. We omit the details. □ 

Remark 14. The homogeneity of the non-local term (in the double-integral form) 
will often simplify our calculations. In this respect, the local nonlinearity may be 
harder to treat than our convolution term. 

The next Lemma takes care of another technical ingredient of the perturbation 
technique. 

Lemma 15. The operator W'^Ieiz^^^): E ^ E is a compact perturbation of the 
identity operator. 

Proof. The second derivative V^/e(z£_^) acts on generic elements v, w G E a,s 

(9) {V'^Ieiz^,^)v,w) = J {Vv -Vw + V{ex)vw) dx - K{v,w), 
where 

(10) K{v,w) — I W{x — y)zs^^{x)^v{x)w{y)dxdy 

W{x - y)ze^^{y)v(y)z^^^{x)w{x) dx dy. 

We need to prove that X is a compact operator. Take any couple of sequences 
{vm}, {wm} in E, such that \\vm\\ < 1, \\wm\\ < 1, Vm vq and Wm wq (weak 



A NOTE ON SCHRODINGER-NEWTON SYSTEMS 



7 



convergence in E). Without loss of generality, we assume vq = Q = wq. We are 
going to show that K{vm, w„i) — > 0. Consider first 



W{x - y)z^^^{y)v{y)zs^^{x)w{x) dx dy 



Here we have chosen p — t in the Hardy-Littlewood-Sobolev inequality. Now, 



(11) Ike.C"™!! 



\z^.^^{x)Vmix) \ ■' dx. 



\Ze,i{x)Vm{x)\'i dx- 
\x-^\<R J\x-i\>R 

Since e H^{B{^,R)), we have Vm strongly in L'^{B{^,R)). The Holder 
inequality implies then 



x-^\<R 



\Zs,^{x)Vmix)\'^ dx < 



Ju 



\Vm{x)\'^ dx 



|2:-5|<_R 



\Ze,i{x)\'^ dx 



o(l) 



as m ^ +0O. To estimate the integral outside the ball B{£_,R), pick 77 > 0: there 
exists i? > such that 



x-i\>R 



\Zs.^^{x)Vmix) \ =5 dx < 



-i\>R 



\Vmix)\^ dx 



\x-i\>R 



2 dx\ <Cr] 



thanks to the exponential decay of z^ ^t and the boundedness of {vm}- The term 
with Wm is perfectly symmetric, and in conclusion 



lim sup 

m — 



W{x - y)zs,^{y)v{y)zs^^{x)w{x) dx dy 
By repeating these arguments, we can show that 



< C?7E 



lim sup 

rn — *+oc 



W{x — y)zg_^{x)'^v{x)w{y) dx dy 
and since 77 > is arbitrary, the proof is complete. 



□ 



Remark 16. The argument contained in the proof just finished, and based on Propo- 
sition [51 will be used several times in the sequel. 

We come now to the most important step of this section: the invertibility of the 
gradient V/e(ze_^) along orthogonal directions to the manifold Z. This will be of 
crucial importance when trying to solve the auxiliary equation ([8]) with respect to 
w. 

The following Lemma is technical, and we recall it here for the sake of complete- 
ness. Its proof can be found in [5J, Claim 1 inside the proof of Lemma 5. 

Lemma 17. There exist R £ {\/s, y/e) such that 



i?<|a:-^|<i?+l 



for euer?/i;_Lspan|ze,e,-g|f 

Lemma 18. For some C" > and any e ^ 1, P\I I^{z^,^) is uniformly invertible 
for every £, E M.^ such that \e^\ < 1 and \\ (FV/e(ze,5))"^ II < C". 
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Proof. By Lemma 1151 we need to prove that in a suitable neighborhood of zero 
there faU no eigenvalues of PV/e(ze^^). The proof is rather technical, and we will 
write down only those steps that require significant modifications to the arguments 
of [H Lemma 5]. To begin with, we recall that (see ([9]) and (fTO)) ) 

{\/'^Ie{ze.^)v,w) = J \/v -X/w + V{ex)vw 

-2 [ W{x ~ y)ze^(^{x)w{x)z^^^{y)v{y)dxdy 

W{x - y)ze^(^{xYw{y)v{y) dx dy. 

Choosing v = w = z^^^ and recalling that Ze,^ solves (O, we get 

{V'l,{z,,^)v,w)= J [V{ex)-V{eO]zl^ 

-2 f W{x~y)z,^^ixfz,^^{yfdxdy<^C\\z,^^f<0. 



Let us abbreviate 

dzi.,^ dzg^^ 9zg,g 



X = span<^ Ze,^, 



Again, by the exponential decay of (see Theorem |9]), X C H^{M.'^). The rest 
of the proof consists in showing the inequality 

(12) {V^I^{z,^^)v,v) >C\\vf foraUu±X. 

We will assume that = 1. Next, for R as in Lemma [T71 we select a cut-off 
function : -> K such that Xi? = 1 in B{^, i?), Xi? = in \ B{^,R+1), and 
|Vxij| < 2. We decompose v = vi + V2, where vi = xrv and W2 = (1 — Xr)'^- ^-S 
in [4] , we can easily check that Lemma [17] implies 

{VI,V2)=0{V^) 

\\vi\\'^ f \Vv,\' + V{ex)\vi\' + OiV^) 

\\V2\\''^ I \^V2? + V{ex)\v2? + 0{^e). 

JVl^\B(^,R) 

To prove p^ . we exploit the bilinearity of the second derivative, and split it as 
follows: 

{V^I,{Z,^^)V,V) = {V'^Ie{Ze,{}VuVi) + {V^ Ie{Ze^i)v2,V2) 

+ 2{V^I,{Ze,^)vi,V2). 

As in the proof of Lemma I15[ we can estimate 

(13) \{V^h{z,,^)vi,V2)\ < \{V1,V2)\ 



c 



W{x - y)zex{x)v{x)ze,^{y)v{y) dx dy 

■(C,i?+1)\B(5,_R) 
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The "far-away" term (V^/e(ze_^)w2, ^2) is easier to study than in ,4J, due to the 
homogeneity of the convolution term. Indeed, 

(14) \(V^I,iz,,^)v2,V2)\^ f \Vv2\^ + V{ex)\v2\'' 

Jk3\s(5,_r) 



2 / W{x~y)zs^^{x)v2{x)z^^^{y)v2{y)dxdy 
W{x - y)z^^^{xYv2{x)v2{y) dx dy 



> \Vv2\' + V{ex)\v2\'-C{R){ \Vv2\'] 

\Vv2\^ + Viex)\v2\'' + 0,{\\v2f) 

>l\\v2r+oiv^). 

Again, we have used the Hardy-Littlewood-Sobolev inequality and the properties 
of The last term to estimate is the one "concentrated" on B{£^,R). Now, 



(V'/e(2.,?)t'i,fi) - Ikill -y Wix-y)z,,^{xyviiyydxdy 
-2 W{x - y)zs^^{x)vi{x)zs^^{y)vi{y)dx dy 



|V«ip + V{eO\vi\' + / [Viex) - V{eO] \vi\^ 

B{^,R+l) JB{i,R+l) 



2 / W{x - y)zs^^{x)vi{x)zs^^{y)vi{y)dxdy 

W{x ~ y)zs^^{xYvi{yY dxdy. 



By the exponential decay of z^^^ and the Hardy-Littlewood-Sobolev inequality, it 
is easy to show that 



(V^4(Ze,c)«l,^^l) = j |V«i|^ + m)l«ir +Oe(l) 

-2 I W{x - y)z^^^{x)vi{x)ze^(^{y)vi{y)dxdy 

W{x - y)ze,^{xfvi{yf dxdy. 



It is very important, at this stage, to highlight that the convolution terms play no 
role, since they do not contain any potential function. This is the main difference 
with respect to |;4 , where the nonlinearity has a potential function K in front. 

The last estimate tells us that, up to an error of size ^/e^ the estimate of 
(V^/£(ze_j)ui, wi) follows from the estimate of the quadratic form 

vi^ j |Vvi|2 + l/(£e)kiP 

-2 I W{x - y)ze^^{x)vi{x)ze^^{y)vi{y)dxdy 

W{x - y)zi,^^{xYvi{yY dxdy. 
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with constant coefficient V{eS^). But this is the second derivative (along the di- 
rection vi) of the Euler-Lagrange functional associated to the limiting equation 
([7]). We can repeat the same arguments of Lemma 5], using of course the non- 
degeneracy for the linearized limiting problem given by Theorem H (and RemarkfTO]) 
instead of the non-degeneracy for the linearized Schrodinger equation, and conclude 
that 

(15) {W^I,{z,,^)vi,v^) > C\\v4^ + o,{l). 
Collecting ([12]), (HJ) and ((15]), and recaUing that 

IkilP + lk2|P = 1 + 0(Vi) = llt'll' + 0(^/i), 
we prove (fT2|) . □ 

Remark 19. The non-degeneracy of the limiting problem is used in an essential 
way to show (fTSl) . 

The auxiliary equation ((5]) is formally equivalent to a fixed-point equation for 
the operator 

(16) S,(w)=w-{PVH,{z,^^)y\PVh{z,,^+w)). 

To find w such that w — Se{w), we will construct a function set on which is a 
contractive map. 

3.1. The contraction principle for S^- As in it is convenient to restrict the 
operator to a proper subspace of the orthogonal space (Tz^ ,= 2^ , so that we can 
get precise estimates on its fixed points. We therefore need to construct a suitable 
subspace of (T^^ ^Z^ that is invariant under S^, and on which S'^ is contractive. 
The idea is to begin with radial solutions of the problem 

r -A?i + m|a;|-"u = /(|x|), |x| > i? 

(17) <^ u(x) = l, \x\=R 

y u{x) — > 0, |x| +00, 

where i? > 0, m > 0, a € (0,2], and /: (i?, -|-cx)) ^ R is positive, and decays 
at infinity. Roughly speaking, the reason why ((TTj) is related to our fixed-point 
argument is that the qualitative behavior of the decaying potential V is essentially 
that of m|a;|~", as described precisely by Lemma [T^ Looking for radially symmetric 
solution of (|17p . and performing the change of variable v{r) — u{r)r, we reduce to 
the problem 

{— u"(r) + mr^°'v{r) — f{r)r [r > R) 
v{R) =R lr = R) 

v{r)r^^ — > (r ^ +oo). 

We highlight that this is precisely equation (21) of (4J, where a term has disappeared 
since in our work n — 3. The homogeneous equation corresponding to (jl8p has a 
general solution spanned by two generators vi and V2- Hence, the general solution 
of ([TT]) is spanned by ui{r) — r~^vi{r) and U2(r) = r~^V2{r). 

For the reader's sake, we recall some more properties of vi and U2, which we 
borrow from U; we recall that we are working in M"^, so that some formulas become 
more explicit. In case a < 2, vi and V2 are modified Bessel functions. If we denote 
by Bj the modified Bessel function of first order, and by Bf the modified Bessel 
function of second order, so that they both solve the ordinary differential equation 

r^y"{r) + ry'{r) - (r^ + i^)y(r) = 0, 



A NOTE ON SCHRODINGER-NEWTON SYSTEMS 



11 



then 

2\/m „ 1 



2-a 



Hence, standard asymptotic properties of the Bessel functions imply that 

2Jrn 

2-a 



vi[r) w r ' exp — r * 



V2[r) « r ' exp r 



As a consequence, 



If a = 2, then 



vi(r) = r 2 



and therefore 



Mi(r) 

M2(?') 



-1- Vl+4r. 

r 2 



r 2 



We recall, for future reference, the following result, see [4]. 

Lemma 20. Let ui, U2 be defined as above, and let Cp be a solution of |i7[ j, where 
f: (i?, +oo) — !■ R is a positive, continuous function satisfying the integrability con- 
dition 

-f C30 



r f[r)u2{r) dr < +oo. 
Then there exists a constant j{R) > such that (f(r) < "f{R)ui{r) for all r > R. 



For any R > fixed, we define the set yVj(i?) as the set of all functions w € E 
such that 

'7(i?)Viui(|x|), if|a;|>i? 
Vi, if < R, 



where ui is the function we have introduced above, and j{R) is the constant of 
Lemma [ini Finally, we set 

r,(i?) ^{weE\\\w\\E< coe, w e w,{R) n (t,^,,^)^} . 

We will choose a good value for the constant cq > later (see the proof of Lemma 



Remark 21. It is possible to select the constants in the truncation so that 



\Ze,^{x) + w{x)\ < 



c 



il + \ex\ 

for all a; e R3 ^nd all w e r,{R). 
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By the last remark, we are authorized to write 
1, 



■w\ 



|2 
E 



W{x " 2/)|ze,^(x) + w(a;)p|ze,^(y) + w{y)\'^ dx dy. 



4, 

Since w _L T^^^Z, it is known that any critical point of /g constrained to the 
manifold {z^.^ + w | it; G r£(i?)} is a solution of ([1]). See for a proof of this 

general fact. 

We fix a value p > such that (see Lemma [T2|) 

^(ea; + e^) > to|x|~" for \x\ > p. 

Lemma 22. For cq large enough, and e sufficiently small, we have 

\\Se{w)\\ < coe, 

for all w G (p) . In addition, is a contraction. 

Proof. Let C be the constant given in Lemma [T51 Moreover, 

for some C" > 0. Choose cq = 2C'C. For w G ^e{p), we have 

VS.iw): v^v- {PV''l,{z,,^)y\PV^I,{z,,^^ + w){v)) . 
We apply PV'^Ie{ze,^) and find 

\\PV^h{z,,^){{VS,(w),v))\\ - ||PV2/,(z,,4)(«) -PV2/,(z,,4 
But, for any wi, W2 G T^{p), 

W2) - D'^Ie{z^^^ +w)iwi,W2) = 

-2 W{x ~ y)zi.,^{x)wi{x)z^,^{y)w2{y)dxdy 



W{x - y)ze^^{xfwi{y)w2{y) dxdy 
+ 2 I W{x-y){z^^^{x)+w{x))wi{x) (z^^^iy) + w{y)) W2iy) dx dy 
+ J W{x-y) {ze,^{x) + w{x))^ wi{y)w2{y) dxdy 

= W{x — y)w{x)wi{x)w{y)w2{y)dxdy 
+ / W{x-y) (^{z^,^{x) +w{x)f - z^,^{xf^wi{y)w2{y)dxdy. 



The first term is estimated by 



W{x — y)w{x)wi (x)w{y)w2 (y) dx dy 



< ci||u;i||i6||w2|lL6||w;| 



E ' 



with q > 8/3. Indeed, we apply the usual Hardy-Littlewood-Sobolev inequality, 
and we get 

< \\W ■ Wi\\^6 \\w ■ W2||^|. 



J W{x — y)w{x)wi{x)w{y)w2{y)dxdy 
Now, 
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For any q > 1, let q' = q/ {q — 1). Then 



1/9 / \ 1/9' 



Now, J |w| 49 is finite, since w € ^e{p)- Choose now q > 8/3 (so that r = 3g/4 > 2); 
there exists a constant C2 such that < C2V{ex). Hence 

\ 1/9 / p \ 1/9 

\w\iA < C3U\w\'\wr^\ 

< C4 \w\^Viex)^ <C5||7«|p/«. 
Similar arguments can be used to estimate the term 

- y) {{ze,£,{x) + w{x)f ~ z^,^{xY'^ Wi{y)w2{y) dxdy 

W{x - y) {w{x)(z^^^{x) + w{x))) wi{y)w2{y) dx dy 



Therefore, 

\\PV'l,{z,MVS,{w),v))\\ < C6||u;||^||«||, 
for some S < 1. If wi, G ^e{p)i then 

\\S,{wi) - S,{W2)\\ < II (PV24(Ze,c))"'||x 

X \\PV^I,{Z,^^) {S,{wi) - Se{w2))\\ 

and by the Intermediate Valued Theorem we find 



||S'e(u;i) - Ss{w2)\\ < C7 l^^max^ ||u;2 + s{wi - W2)||J ||wi - W2||- 

Since wi and W2 belong to T^{p), we easily deduce that 

(19) \\S,{wi) ~ S,{W2)\\ = 0,{1)\\WI - W2\\. 

Hence is a contraction on r^{p). Taking wi = w and = in we get 

||5e(u;) - 5,(0)11 -o,(l)||,i;||. 

On the other hand, 

115,(0)11 = II (PV24(zs,c))"'(PV4(ze.c)) 
<C'||FV/,(ze,e)|| <C'Ce. 
By the triangular inequality, 

||5,(«;)|| < \\S,{w) - 5e(0)|| + ||5,(0)|| < 0,(1)1^11 + C'Ce 
= Oe(l)||w|| + ^coe. 

Since w G F,, \\w\\ < cqe; hence ||5e(w)|| < cqE provided e is small enough. □ 

The most important, and probably the most difhcult step, is to prove that 
maps ^e{p) into yV^{p). Our proof follows [1], but we remark that the form of the 
convolution term can be exploited to cancel several terms. 

Lemma 23. For all e small enough, Se(rs{p)) C yVe(p). 
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Proof. For simplicity, wc set w = Se{w). Wc notice that 

W = W- {PV^IeiZe,^))'^ Pyie{Ze,i + w), 

which we can write as 

P {V^Ie{z,,^){w -w)+ V/e(Ze,J + W))=0. 

Hence w satisfies the identity 

— Aw + V{ex)'w ~ ''^ J W{x — y)ze,c^{y)^ dy 

-2zs,^ j W{x-y)ze,^{y)w{y)dy 
= '^w j W{x-y)zs^^{y)w{y)dy + {zs,i+w) j W {x - y)w{yf dy 

- (^Aze,i + V{ex)ze,i- Ze,i j W {x - y)z,^^{yf dy^ 

+ (-A4,5 + V{e)ze,i) V, 

where 

77 = (V2/e(Ze,4)(w -W)+ VIe{Ze,i + w),Ze,^) ps.dl"^ 

and ie^j stands for a linear combination of the derivatives 

d d d 

related to the projection of the equation w = Se{w) onto (Tz^ a^)^ ■ 

We define zo{x) = Zs,^{x + $,), wo{x) = 'w{x + $,), and ^0(2;) = w{x + Let Lq 
be the linear operator defined by 

Lqv = —Av + V{ex + e^)v ~ J — y)zo{y)^ dy 

- 2zo j W{x - y)zo{y)v{y) dy. 

Therefore wq verifies 

LqWo = 

2wo j W{x-y)zo{y)wo{y)dy + {zo + wo) j W {x - y)wo{yf dy 

- i-AzQ + V{ex + e$)zo - zq W{x - y)za{yf dy 



+ {-Azo + V{ex + eC)ZQ) T]. 

If we set 

gi = 2wo j W{x-y)zQ{y)wQ{y)dy+{zo + wo) j W{x - y)wo{y f dy 
g2 = {-Azo + V{ex + e$,)zo) r? 

53 = -Azo + V{sx + s^)zo - ^0 J W{x - y)zo{yf dy, 

and go = g\ + g2 + 93, then Lowo = go- 
We need to prove that 

(20) \wo{x)\ < \fe for \x\ < p 
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and 

(21) \wo{x)\ < ^{p)V^ui{\x\) for |a;| > p. 
The main ingredients are of the foUowing estimates: 

(22) I15o||l«/2(b(o,2p) < Oe(l)V£ for <? > 3 
and 

(23) |5o(a;)| < o,{l)y^ui{\x\) for \x\ > p. 

Moreoever, |gi(a;)| < ci|wo(a;)|. Indeed, this follows easily from the definition of gi 
and from the inequality 



Wi- -y)zo{y)wQ{y)dy 



< C2 1 1 2:0^0 1 1^3/2, 



which is the limit case of ([5]) as s ^ +00. Now (22) and follow from arguments 
that are completely analogous to those carried out in 4J, pages 338 and 339. We 
just suggest how ^ and jH]) are deduced from ^ and ((23l) . 

By standard elliptic estimates (see for instance [13], Theorem 8.24) we have, for 

g>3, 

||wo||l=°(B(0,p)) < C3||wo||l2(B(0,2p)) + C4||5o||l9/2(b(0,2p)) • 

Since ||2:e,{|| < C5 and by virtue of Lemma we deduce 

\\wa\\L-^B{o,2p)) < ces. 
Coupling this with (22), we find 

Hence (PUI) follows. If a < 2, the proof of (pij) follows as in by comparing wq 
with the solution of the problem 

-A(j}+\x\-"<j)^ ^/eui, \x\>p 
(t){x) = = P, 

exploiting Lemma 1201 A few changes are in order when a = 2, and we refer to 

m- □ 

We have all the ingredients to state a result about the existence and uniqueness 
of a fixed point for Se in T^^p). We omit the proof, which is based on standard 
techniques and can be found in O HI [5] . 

Proposition 24. Under the assumptions of the Theorem^ for s sufficiently small, 
there exists one and only one fixed point w € Fg (p) of ■ Moreover, w is differen- 
tiable with respect to S M^, and there exists a constant C > such that 

(24) ||w|| < Ce, llV^wll < Ce\ 

where 6 < 1 was found in the proof of Lemma\2^ 

4. Analysis of the reduced functional 
We recall that the manifold 

^ = + w I e w e Te{p)} 

is a natural constraint for the functional : this is a general feature of the perturba- 
tion method that we are using, see |3j. Loosely speaking, we have found w _L T^^ 
such that V/e(ze,5 + w) G ^Z. If /g has a critical point constrained to Z, then for 

£ sufficiently small we have V4 € T^^^^Z r^ {T^^^Z)^ = {0}. In addition, Z ~ , 
and therefore the constrained function given by 

$e(0 -4(^e,C + U;) 
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has the property that every critical point ^ produces a solution z^^^ + w of (|3]) for 
£ sufficiently small. However, a direct study of <i>e on R'^ is difficult, and we want 
to use its Taylor expansion to locate its critical points. 

Proof of Theorem^ We expand the function <i>e: 

'^eiO = \j\V{z,,^+w)\^ + \jv{ex)\z,,^+w\^ 

If 2 2 

- ^ I W{x-y){ze^^{x)'Vw{x)) {z^,^{y) + w{y)) dxdy 

^ ^ ' — |2 I Iv7„..|2 I I v7„ V7„., I ^ I ^rt — 



2 / ' 2' 



^ 7 / ^ i^eAx) + w{x)f {ze,^{y) + w{y)f dx dy 



W{x - y) (ze,j (x) + w{x)f {ze,i{y) + w{y)f dx dy. 



Since 



lk..dl'+ / iVi^O-V{£x))zl^^ f W{x-y)z,,^ixfz,,^iyfdxdy, 



{z^.^,w)e + J (Vie^) - V{ex)) z^^^w 

= j W{x - y)z^^^{xfze^(^{y)w{y)dxdy 

we have 

j W{x-y)z,.^^{xfz,,^{yfdxdy + ]^ j {V {e^,) - V {ex)) zl^ 
+ j W{x-y)z^^^{x)'^Ze^(^{y)'w{y)dxdy+'^\\w\\'^ 



/ W{x-y){ze^(^{x)+w{x))'^ {ze,^{y)+w{y))^ dxdy 



+ j {V{eO-V{ex))z,,^w 
= i / W{x - y)z,^^{x)^z,^(^{y)^ dxdy + 'Re{w), 



where 



T^e{w) = - / {V{ei)-V{ex))zl^ 



+ W{x - y)z^^^{x)^Ze^^{y)w{y)dxdy 



If 2 2 

W{x-y){zs^^{x)+w{x)) {z^,^{y) + w{y)) dxdy 

+ / {V{eO~V{ex))z,^^w + ^\\w\\'=o{e) 

by the usual arguments based on the Intermediate Valued Theorem and on the 
Hardy-Littlewood-Sobolev inequality and by Proposition [M] The calculations are 
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completely algebraic, and a direct use of the estimates proves the statement. 
We omit the routine calculations. Inserting the definition 

z,^^{x) = V{eC)U (ym) (x - O) , 

into the expansion just obtained gives 

(25) ^,{^)^^V{eO^ j W{x-y)U{xfU{yfdxdy + o{e). 

If (we may of course assume that xq = 0) is an isolated local minimum (or 
maximum) of V , then ^e{C) — possesses a critical point ~ 0, and this 

gives rise to a solution 

...(f)--..,(^) 

concentrating aX Q — xq. □ 

Remark 25. Since z^^^ is positive, and w — w{e,^) is of order e, the solution Ug 
found in Theorem |4] is positive. 

5. Other existence results and comments 

Although wc have proved our main theorem under the assumption that V has an 
isolated minimum or maximum point, some extensions can be proved. Since their 
proofs require some rather standard modifications (the statement about a compact 
manifold of non-degenerate critical points was studied in [S], see also [31ll]) of the 
perturbation technique, we collect these generalized existence results in the next 
theorem. We recall that a stable critical point xq of a smooth function is a critical 
point at which the Leray-Schauder index ind(Vy, xo,0) ^ 0. It is well-known that 
strict local minimum (and maximum) points are stable critical points. 

Theorem 26. Theorem is also valid provided xq in an isolated stable critical 
point of V. More generally, if S C M"^ is either a compact set of maxima (or 
minima) of V , or a compact manifold of non-degenerate critical points of V\^^± , 
then the conclusion of Theorem^is again true. 

As recalled in the first section, in [12| a more involved system with a magnetic 
field has been recently studied. It has been known since the paper [11] that the 
electromagnetic Schrodinger equation can be handled by means of the perturbation 
method discussed in our paper. However, it generally requires some boundedness 
assumptions on the potential that generates the magnetic field, and this excludes 
many fields of physical interest. However, we believe that our Theorem S] holds 
(and we find solutions such that \ue\ concentrates at "good" critical points of V) 
for the non-local equation 

(^jV - A{x)y u + V{x)u = (VF* \u\^)u 

if the following technical assumption is satisfied (along with (VI) and (V2)): 

(Al): A: R'^ R'^ is a bounded vector potential, with bounded derivative: 

SUPa;GR3 < OO. 

Of course solutions are complex-valued, but their moduli still concentrate. The 
action of A is felt only in the phase of these solutions. This is, loosely speaking, a 
generic feature, as proved in [20]. To feel the presence of A, one usually needs some 
singularities in the magnetic field, or a very strong symmetry of the solution. See a 
recent paper by Cingolani et al. ([TOj) for the latter case, though in the local case. 
We plan to investigate the non-local case in a future paper. 
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An interesting open (as far as we know) problem is the generalization of our 
results to any dimension rt > 3 (the convolution kernel W becomes homogeneous 
of degree 2 — n, of course). Even for a potential V bounded away from its infimum, 
there are obstacles. Every paper in the literature deals with n — 3. We cannot say 
whether this is only a technical obstruction. With respect to the method we have 
used, the main difficulty is to extend Theorem [HI in particular the non-degeneracy 
part. We suspect that the ODE proof contained in [22^ can be adapted to any space 
dimension. 

Another interesting question is whether potentials V with a fast decay can be al- 
lowed. More precisely, this corresponds to the assumption that liminf i^-i^^oo V^(a;)|xp = 
0. The only result which deals with such potentials is [TH], but, once more, it seems 
difficult to carry their penalization scheme to our non-local situation. Moreover, as 
stated in [TH], it would be useless to look for solutions near our manifold Z^, since 
one expects solutions to decay polynomially fast at infinity. 
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